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We propose the two-band s-d model to describe theoretically a diffuse regime of the spin-dependent 
electron transport in magnetic tunnel junctions (MTJ's) of the form F/O/F where F's are 3d tran- 
sition metal ferromagnetic layers and O is the insulating spacer. We aim to explain the strong 
interface sensitivity of the tunneling properties of MTJ's and investigate the influence of electron 
scattering at the nonideal interfaces on the degradation of the TMR magnitude. The generalized 
Kubo formalism and the Green's functions method were used to calculate the conductance of the 
system. The vertex corrections to the conductivity were found with the use of "ladder" approx- 
imation combined with the coherent-potential approximation (CPA) that allowed to consider the 
case of strong electron scattering. It is shown that the Ward identity is satisfied in the framework 
of this approximation that provides the necessary condition for a conservation of a tunneling cur- 
rent. Based on the known results of ab-initio calculations of the TMR for ballistic junctions, we 
assume that exchange split quasi-free s-like electrons with the density of states being greater for 
the majority spin sub-band give the main contribution to the TMR effect. We show that, due to 
interfacial inter-band scattering, the TMR can be substantially reduced even down to zero value. 
This is related to the fact that delocalized quasi-free electrons can scatter into the strongly localized 
d sub-band with the density of states at the Fermi energy being larger for minority spins compared 
to majority spins. It is also shown that spin-flip electron scattering on the surface magnons within 
the interface leads to a further decrease of the TMR at finite temperature. 

PACS numbers: 75.70.-i, 73.40.Gk, 73.40.Rw, 85.30.Mn 



I. INTRODUCTION 

Magnetic tunnel junctions (MTJ's) consisting of two 
ferromagnetic layers separated by the insulating spacer 
exhibit the tunneling magnetoresistance effect (TMR) 
when they are switched by applying a magnetic field from 
the antiparallel to parallel alignment of magnetizations 
in magnetic layers. The TMR effect was first observed 
by M. JuUiere in 1975EJ. JuUiere found that the tunnel- 
ing conductance of the trilayer structure Fe/Ge/Co de- 
pended on the angle between magnetizations in the Fe 
and Co layers. The measured amplitude of the TMR 
in these experiments was 14% at 4.2K. Only 20 years 
later, the large values of TMR at room temperature were 
obtainetLjin magnetic junctions comprising the AI2O3 
barrieraEla. Since that time, there has been a renewed 
technological and fundamental interest to the tunneling 
phenomenon and during the last decade a lot of exper- 
imental and theoretical papers were published on this 
topic (see reviews [5,6,7,8]). pgr-, 

Experimentally it was observedB'B'El that the TMR de- 



pends critically on the material of the insulating bar- 
rier and on the conditions of its preparation, in partic- 
ular on the imperfections of t|ie-4nterfaces between the 
metal and the insulating layemlld. On the other hand, 
the first theory suggested by M. JuUiere expresses the 
TMR ratio in terms of the effective spin polarizations 
Pi and P2 of two magnetic layers via the expression 
TMR = 2PiP2/{l + P1P2), and thus predicts no depen- 
dence of the TMR on the parameters of the barrier. In 
spite of JuUiere's formula is widety used for the interpre- 
tation of the experimental dajtaO, it does not have rig- 
orous theoretical foundationtll. The dependence of the 
TMR magnitude on the parameters of the tunnel bar- 
rier and the metal/insulator interface was takeftjinto ac- 
count in the subsequent theories. Slonczewskll3 consid- 
ered a quantum mechanical problem of tunneling of a 
free electron through a simple spatially uniform barrier 
and showed that the TMR ratio depends on the height 
of the potential barrier and on the effective mass of the 
tunneling electron inside the insulator. Later on, the 
influence of spin-flip scattering at the interfaces on tem- 
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perature and bias-voltage dependences of the TMR was 
investigated by Zhang et al.cB. It was shown that mix- 
ing of spin-up and spin-down tunnel channels leads to a 
decrease of the TMR. 

Due to a permanent progress in the development of 
different methods of ab-inito calculations of the elec- 
tronic properties of solids, in the last four years the 
transport properties of tunnel junctions were investigated 
using the realistic band structure of the ferromagnetic 
layers and the insulator (Refs. [14-20]). The systems 
where the conductance—was calculated were ideal epi- 
taxial Fe/2^e/Fe(001)E3 and Fe/MgO/Fe(001) tunnel 
junctionst3E2l. For these structures the electron trans- 
port is assumed to be ballistic, i.e. the electron momen- 
tum fc|| parallel to the ferromagnet/insulator interface is 
conserved. Experimentally, it was also demonstrated re- 
cently by Heinrich et alH that Fe/MgO/Fe(001) junc- 
tions can be indeed grown by depositing the MgO epi- 
taxially onto a Fe whisker and then depositing another 
Fe electrode epitaxially on a top of the MgO. The re- 
sults of ab-initio calculations showed that tunneling has 
more complicated behavior than the predictions, of the 
simple barrier model proposedj-biuSlonczewskill3. The 
main conclusions are as foUowsEaEj. (i) Tunneling con- 
ductance depends strongly on the symmetry of the elec- 
tron states in the ferromagnetic electrodes and in the 
insulating layer, (ii) The decay rates of evanescent states 
in the barrier are different for the states with different 
symmetry. The slowest decay rates have the evanescent 
states which are compatible with s symmetry. The Bloch 
states in the metal couple more efficiently through the 
interface with the decaying states of the same symmetry 
in the barrier. Therefore, mostly free-like sp-electrons 
from the bands with s character in the ferromagnet give 
the essential contribution to the tunneling conductance, 
(iii) For thin insulating layers the tunneling current in 
the minority channel is dominated by the interfacial res- 
onance states that gives rise to "hot spots" for the fc||- 
rcsolved conductance in the surface Brillouin zone. The 
contribution from the resonance states is substantially 
suppressed for the thicker barriers sipce the Bloch states 
at "hot spots" have no s chacacteilla. Moreover, as it 
follows from recent discussion£3, the surface resonance 
states are very sensitive to the details of the interface. 
The asymmetry of the potential barrier and the interfa- 
cial roughness considerably reduce the contribution from 
the surface states to a total conductance. 

Nowadays the ab-initio calculations of the TMR are 
possible only for ballistic junctions like above mentioned 
Fe/ZnSe/Fe and Fe/MgO/Fe structures which are char- 
acterized by a small lattice mismatch between the metal 
and the insulator, and by the well-defined band structure 
of the oxide or the semiconductor. In realistic junctions 
the electron transport has a diffuse character, i.e. when 
the electron crosses the ferromagnet/insulator interface 
its fc|| momentum is not conserved. For F/O/F tunnel 
junctions the mppt successful material until now has been 
alumina, Al203a. The AI2O3 tunnel barriers are usually 



fabricjated by the natural or plasma oxidation of the Al 
layeitfla. The subsequent structural analysis, e.g. with the 
use of the X-ray photoelectron spectroscopy, shows that 
alumina is amorphous and the obtained AIQ— tunnel bar- 
rier deviates from the ideal AI2O3 structuret2l. For a uni- 
form coverage the Al film thickness is usually ranged from 
about 7 tQ 18 A, depending on a type of the ferromagnetic 
electrodeO. There is a small range of Al thicknesses that 
yields to the best TMR ratio for a given oxidation condi- 
tion. When the Al layer is thin, the ferromagnet surface 
becomes oxidiaed leading to the formation of CoO and 
C03O4 oxideso or to the Fe304 oxidetJ. On the other 
hand, with too thick Al film, an excess of Al metal is 
left unoxidized. The amorphous barriers, the roughness 
of the interface and its structural inhomogeneity make 
the rigorous ab-initio calculations of the TMR virtually 
impossible and therefore more simplified models are re- 
quired to treat the diffuse electron transport in MTJ's. 

In our previous papeiS, we attempted to investigate 
the influence of scattering processes of the electrons (with 
and without spin-flip) at the interfaces on the TMR using 
a simple two-band (spin-up and spin-down) free-electron 
model. It was shown that for this simplified model spin- 
conserving scattering may or may not lead to a decrease 
of the TMR depending on the amplitude of the scatter- 
ing potential. In this paper, we proceed to study the 
diffuse electron transport in magnetic tunnel junctions 
of the from F/O/F with F's being 3d transition metal 
electrodes and O being the insulating barrier (AI2O3). 
We use the results of ab-initio calculations to model the 
band structure of ferromagnetic electrodes, namely, the 
most important feature of these calculations that at least 
two groups of electrons form the total band structure: 
the almost free-like spin-up and spin-down bands and 
the narrow strongly exchange split bands. We will call 
the first group s and the second group d electrons and 
will consider that a periodic part of the s-d hybridiza- 
tion between bands is taken into account that results in 
the nonequivalence of spin- up and spin-down s-bands. 
The parameters of the adopted model can be adjusted 
to reproduce a value of the TMR observed in the experi- 
ments. We assume that for the case of AI2O3 barrier the 
exchange split s-like quasi-free electrons give the main 
contribution to the TMR effect. In the framework of this 
simplified model we will be able to investigate in a proper 
way the influence of electron scattering at the interfaces 
on the tunneling conductance. 

To describe the nonideal tunnel junction we assume 
that defects and impurity centers (the Al or O ions, 
or other artificially embedded ions) are randomly dis- 
tributed within few monolayers near the F/0 interface. 
Within the interface, an electron undergoes scattering 
when it comes to the defect or impurity center. We take 
into account these processes assuming that the hybridiza- 
tion between s and d bands changes randomly on the in- 
terface because of this parameter is the most important 
in the adopted model among the other ones, character- 
izing the scattering potential. It yields to the possibility 
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of s-like electrons to scatter into the d sub-band (and 
vice versa) and thus strongly affects tunneling. To treat 
the electron scattering on the random potential we use 
the coherent-potential approximation (CPA)lJ that al- 
lows to consider the case of strong scattering. We apply 
the generalized Kubo formalism and the Green's func- 
tions method to calculate the conductance of the system, 
and find the vertex corrections to the conductivity with 
the use of "ladder" approximationEa combined with the 
CPA. It is shown that the so-called Ward identity is satis- 
fied in the framework of this approximation that provides 
the necessary condition for a conservation of a tunneling 
current. Note, that if it is not the case, the conclusions 
may be completely misleading. As a result, we show that, 
due to substantial difference in majority and minority d 
density of states at the Fermi energy for 3d ferromagnetic 
metals, the inter-band s-d scattering on the interface can 
strongly reduce the TMR even down ta.zero value. In ac- 
cordance with results of Zhang et alS3, it is also shown 
that the spin-fiip scattering of electrons on the surface 
magnons within the interface leads to a further decrease 
of the TMR ratio at finite temperature. 

The paper is organized as follows. In Sec. II we de- 
scribe the model Hamiltonian, the calculation of the tun- 
neling conductance and vertex corrections. The discus- 
sion of the obtained results is presented in Sec. III. Con- 
clusions are in Sec. IV. The proof of Ward identity is 
given in Appendix A. The details of the derivation of the 
CPA equations are described in Appendix B. 



II. THEORETICAL MODEL 

A. The Hamiltonian of the system 

We will consider a trilayer tunnel junction of the form 
F1/O/F3, where Fi and F3 are two semi-infinite fer- 
romagnetic layers and O is a dielectric oxide spacer 
(AI2O3). Our arguments in behalf of the two-band s-d 
model which was briefly described in the introduction are 
as follows. One of the conclusions of ab-iniliojpalculations 
of the TMR for ballistic tunnel junctionstfEj is that the 
expected spin-dependence of the tunneling current can 
be deduced from the symmetry of the Bloch states in the 
ferromagnet at the Fermi energy. Spin-polarized band 
structure for bcc Fe, fee Ni, and fee Co can be found in 
Ref. [26] . The type and symmetry of the Bloch states for 
different crystal faces with k|| = for Fe, Co and Ni are 
presented in Table I (in accordance with Ref. [18]). For 
example, in case of Fe electrodes, the examination of the 
band structure shows that both the majority and minor- 
ity bands with s character in [110] and [111] directions 
(Ai and Ei) cross the Fermi energy. For [100] direc- 
tion the band with s symmetry (Ai) crosses the Fermi 
energy in case of the majority channel only. The simi- 
lar analysis can also be performed in case of Ni and Co. 
Thus one can assume that in the polycrystalline Fe, Co 
or Ni-based films the states with s character present for 



TABLE I: Type and symmetry of the Bloch states with 
k|| = for Fe, Co and Ni for three different crystal faces 
(in accordance with Ref. [18]). The symmetry of these bands 
is as follows: Ai, Ei, and Ai {s,p,d); A5, and E2 (p and d); 
and A2, A2', E4, and A3 (d). 



100 



110 



111 



Fet 
Fei 
Cot 
Coi 
Nit 
Nil 



Ai, A2', As 
A2, A2', As 

Ai 
Ai, As 

Ai 

Ai, A2, As 



El, E3 
El, E3 

El 
E2, E4 
El, E3 
El, E2 



Ai 
Ai 



both spin directions and the Bloch states will couple ef- 
ficiently through the F/0 interface with s states in the 
insulator and will decay with the equal rates in the bar- 
rier region. We will call the electrons from these bands 
as s-like electrons and will describe them as free elec- 
trons with the effective mass rris ~ trie (where me is a 
bare electron mass) and with different Fermi momenta 
kf^ and kf^ {kf^ > kf^) for up and down spins. The 
idea about a dominant contribution of the mostly itin- 
erant electrons to tunneling was originally proposed by 
M. Stearn£0 and explained the positive polarization of 
the spin-dependent current in tiie experiments on tun- 
neling with the superconductoraJ. According to estima- 
tions of StearnsO, in the case of Fe, kP = 1.09 A-i, 
kf^ = 0.42 k-\ 

Other more localized bands (compatible with d sym- 
metry) crossing the Fermi energy also will be described 
by two exchange split bands with the isotropic quadratic 
dispersion law but with larger effective mass ^ m^. 

FT F I 

The Fermi momenta k^ and k^ of d-like electrons can 
be chosen to reproduce the typical for 3d transition met- 
als ratio of the values of spin-up and spin-down d density 



of^ates at the Fermi energy £i?, P^(£f) : Pdi^p) ~ 1 : 



H-accordance with the band structure of Fe, Co 
and Nit£IE3 the narrow majority and minority c?-bands are 
practically filled. Therefore, the particles with a large ef- 
fective mass, TTid, must be regarded as holes. The values 
of the Fermi momenta define the positions of band bot- 
toms (Fig. 1) 1/-^, y.^ [a = s,d; i — 1,3) with respect 
to the Fermi energy ep- We note, that the aim of this 
work is to calculate the relative change of the TMR due 
to scattering. Therefore the proposed model is rather 
adequate for this purpose since the scattering rate de- 
pends mostly on the density of states and on the matrix 
elements of the scattering potential. 

The AI2O3 tunnel bajpders obtained by oxidation of the 
AI film are amorphoustltj. Concerning the a-Al203 cryjSa 
tals, it is known from the band structure calculationsE3 
that the gap (which is not direct) between the upper va- 
lence band and the conduction band is of the width of 
w 6.29 eV. The dispersion law of the lowest conduction 
band is not isotropic, and the effective electron masses 
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FIG. 1: The model potentials describing the propagation of 
an electron in the trilayer tunnel junction F/O/F. The solid 
Une corresponds to the potential profile of s-like electrons and 
the dashed line for d-like electrons. V^^^^ denote the spin- 
dependent band bottoms, Ua and Ud are the bottom and the 
top of the conduction and valence bands in the insulator, ef 
is the Fermi level, zo is a typical width of the interface (detail 
description is given in the text). 



along the different directions in the Brillouin zone vary 
from O.IB rue to 0.40 me with an average value of about 
0.35 mec3. This lowest conduction band is formed by a 
mixture of Al 3s, O 2s and O 3s orbitals. In the top of the 
upper valence band the dispersion curves are very flat, i.e. 
the effective masses of holes are large as compared with 
the mass of the conduction electrons. At the F point, the 
top of the valence band consists of hybridized O 2p and 
Al 3p orbitals. In a view of that, one can describe the 
conduction and valence bands of the amorphous AI2O3 
by isotropic quadratic laws with effective masses and 



1771° I »m°. 



Next, the s-like electrons from the ferromagnet 
can easily penetrate into the oxide since the lead 
bands with s character can couple efficiently with 
s states in the barrier, and the tunneling conduc- 
tance, caused by the specular transmission of the Bloch 
waves through the interface, will decay for both spin- 
up and spin-down s-channels with the same law as 



exp 



'2w 



where w is the 



width of a barrier. Us is a bottom of the conduction band 
and K = fc|| is the electron momentum xiarallel to the F/0 
interface. Following MacLaren et alS^ we can suppose 
that d-like electrons from those bands in the ferromag- 
net without s character can not couple efficiently with s 
states in the oxide. Nevertheless, these d-bands have an 
admixture of p symmetry and therefore can couple with 
the valence bands of the AI2O3. However, in the d chan- 
nel the tunneling conductance due to specular transmis- 
sion will decay very rapidly with the thickness w of the 

,^{2ml/h^)ieF'Ud) + K^\ , be- 



barrier, as ~ exp 



cause of a large mass of holes in the alumina (here Ud 
is the top of the upper valence band). The model poten- 



tials describing the propagation of the electron through 
a tunnel junction are presented in Fig. 1. 

We suppose that the ferromagnet/insulator interfaces 
are rough and contaminated by impurity ions. The con- 
tamination of a few metal monolayers close to the F/0 
interface always presents after thej-cpcidation ofpthe Al 
film. Other impurities (e.g. Cr, RuE3 or Cu ion£J) can 
be artificially inserted into the interface. We will charac- 
terize these structural defects by the random potentials, 
which may be divided into the spin-conserving and the 
spin-flip parts. Let zi — a and Z2 — b he the positions 
of interfaces, the xy-plane be the plane parallel to the 
interface, and the z-axis be the axis perpendicular to the 
barrier (see Fig. 1). The Hamiltonian of the system is 
written as: 



H Hq -f^spin— cons. ^" -^^spiii— flip ; 



(1) 



with 



Ho =J2 [ drr^+{r) 
+ J2 J rfrVf (r) 



2tos(z) 



2md{z) 



V^ + WAz) 



V' + U'Jz) 



spin— cons. 



E EE / dr5{z-Ze,)S{p-pZK 



-^spin— flip 



f dv5{z-z^)5{p-p:,)j: 

X [^^+{v)^{v)s.{pf^+r^{v)r^{v)s+{pi) 

Here =t , i is a spin suffix; ipp'{r), ^J^+(r) and V'^(r), 
-0^ (r) are field operators of the creation and annihilation 
of s and d-type electrons with the spin p at the point r; 
ms{z) and md{z) are the effective masses of electrons 
and holes in the ferromagnetic layers (m^ and rud) or 
in the barrier (m^ and m^). U^^{z), Ul^{z) are the spin- 
dependent step- like potentials shown in Fig. 1: 




z < a + Zq 
Ui{z) ^{ Ud a + zo<z<b- 
z > b - Zo, 



Zo 



(2) 



^3d 



where Vf^, 



1,3) are the majority {]) and mi- 
nority (I) band bottoms in the ferromagnet, U^ and Ud 
are the bottom of the conduction band and the top of 
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the upper valence band in the insulator. The positions 
of "steps" for the potential U'^{z) are shifted by a value 
zq with respect to the points z = a and z — h in order to 
describe the finite thickness (~ zo) of the interface. The 
explanation, why it is done in this way, is given below in 
the text. 

^?^spin-cons. is the spin-conserviug part of the Hamilto- 
nian, a = 1, 2 are the interface numbers. To describe the 
defect structure of the nonideal F/0 interface we consider 
that the impurity ions and the ions of a ferromagnetic 
metal constitute the more or less random alloy of a type 
Aa;Bi_2; where A denotes the ions of the ferromagnet (Fe, 
Co or Ni) and B denotes the impurities. Therefore, we 
suppose that each site on the interface a is charac- 
terized by the random parameter 7" of s-d hybridization 
taking two different values and 7s with the probabil- 
ities X and (1 — x), respectively. We also suppose that a 
periodic part of s-d hybridization is taken into account 
leading to the nonequivalencc of majority and minority 
s-bands in the ferromagnet. 

For the sipaple two-band tight-binding model of the 
binary alloysO one assumes that: 1) parameters Eg and 
Ed describing the positions of s and d bands on the energy 
scale are different for the alloy's components; 2) also the 
parameter 7"^^^) of the hybridization between s and d 
bands depends on the type of an ion (A or B). In our 
particular case, one can not take into account the former 
effect since the well-defined two types of the electrons 
will as before exist in the vicinity of the interface and the 
adopted model of tunneling will not change significantly. 
However, the random hybridization makes possible the 
processes of scattering of quasi-free s-electrons into the 
localized d sub-band and vice versa and, therefore, can 
strongly influence on tunneling. 

^^spin-flip is a part of the Hamiltonian describing the 
spin-flip scattering. We take into account only the spin- 
flip processes for s-like electrons since these electrons are 
itinerant and give the most essential contribution to the 
tunneling current. Operators 5^ (p"), S'~(p") are defined 
as 

q 

Here 6q and 6q denote the creation and annihilation op- 
erators of the surface magnons, iV is a number of the 
lattice sites on the interface, 5" is a spin value. We used 
the well-known representation of the spin operators in 
terms of 6q and 6+ and left the first terms of the series. 
J" is a random exchange integral which also takes the 
values J A and Jb with probabilities x and {1 — x). 

Let us now turn back to the step-like potential for d- 
holes. As it will be clear from the following consideration, 
the amplitude of the effective scattering potential on the 
interface (z = a) for s-like electrons (i.e. — Iml]*^(a), 
where S"* is the self-energy) is determined by the value 



a;(l - x){"fA - 7B)^IniG''''(a), here Im G"^''(a) is the 

imaginary part of the c?-electron Green's function at the 
point z — a (density of states). If one put zq = in 
Eq. (H), then the vertex contribution to the conductiv- 
ity (that is the contribution due to tunneling assisted by 
interfacial scattering) will be rather small as compared 
with the contribution due to direct tunneling. This result 
is not enough accurate and is the sequence of a contin- 
ual type of the model when one neglects the existence of 
the atomic lattice and for this reason the self-energy has 
a (5-like behavior on the interface. However, the width 
zo of the interfacial layer is about the distance between 
atomic planes or even larger. For the case of bcc lattice, 
^0 = ao/2 for [100] direction (here ag is a lattice con- 
stant). One can show that for the present model with 
Zo = the imaginary part — ImGQ''(z = z') of the un- 
perturbed d-electron Green's function (with I]'*''(a) = 0), 
i.e. density of states of d-electrons (holes), has small value 
at the point of interface z = a and increases inside the 
F-layer up to the distance of the order of zq ^ ag, and 
then it oscillates near the average value which is approx- 
imately ten times larger than 



ImGg'^(a). The period of 



oscillations is determined by fc^ or fc^ , depending on 
the electron spin and on the orientation of magnetization 
in the F-layer. Such a behavior is easily understood, if to 
take into account that d-electrons (holes) are almost com- 
pletely reflected on the metal/insulator interface. Thus, 
if one shifts the positions of "steps" for U^{z) with re- 
spect to U^{z) as it is given by (y), one can expect the 
more effective mechanism of scattering due to the larger 
value of — Im G^"^, and therefore the model becomes more 
realistic. 



B. The calculation of the tunneling conductance 

To calculate the non-local conduotivity we apply the 
Kubo formula of the linear responses (it is valid under 
the small applied voltage which is much less than the 
value of the energy gap in the insulator): 



a^p(r,r') = ^ Sp { (G;?,(r, r')^r^r'G^^(r', r) 
+ (G;tp(r,r')i)r^r'G^^(r',r))}, 

where a matrix operator is defined as 

/ T7_ 

1 



(3) 



2ms (z) 




(4) 



2md(z) 



and Vr== (Vr — Vr) is the asymmetric gradient operator, 
G^p(r,r') and G^p{r,r') are the retarded and advanced 
(2 X 2)-matrix Green's functions (with components ss, 
sd, ds and dd), =t)i are the spin suffixes, brackets 
(. . .) denote the averaging over the configurations and 
magnon degrees of freedom, the trace (Sp) goes over s 



6 



and d indices of the bands. Below, for convenience, it is 
assumed that Ti = 1. We will recall about h in the final 
expressions for the conductance. 

To calculate the conductivity of the system one has 
to find the Green's function of the Hamiltonian (|l|) , which 
can be found by solving the following system of |differen- 
tial equations in the mixed (k, z) representational: 



E 



[e^"7 _ H-^z)] Gli{z,z') = 5^p5{z ~ z'\ (5) 



along z-axis in Fi and F3 layers, respectively (k is the in- 
plane component of the momentum, a is a band index), 
and let = \/ (^"q)^ + be the imaginary electron 
momentum in the insulating layer. 

Further, for convenience, we will omit indices /i and k 
in the notation of some functions. Equation has four 
linear-independent solutions which we denote as 




(z-1,2) 



Hl^z) 



+ + '^"^ (6) 



1 

E^^(a)<5(z-a) + - 6), 



and 



where k — fc|| is a projection of the electron momentum 
on the xy-plane (parallel to the F/0 interface) , the Greek 
indices a, /3 and 7 assume values s and d, [i denotes the 
electron spin. H^^{z) is the (2 x 2)-matrix linear dif- 
ferential operator where S^^(a) and T,'^^{b) {a, (3 = s,d; 
M =Tii) denote the coherent potentials for spin-up and 
spin-down electrons which take into account the scatter- 
ing of the electron by random spin-conserving and spin- 
flip potentials on the interfaces. They were found witli 
the use of the coherent potential approximation (CPA)EJ, 
the details of these calculations are presented in the sub- 
sequent section IIC and in Appendix B. The operator 
H^{z) represents the effective single-particle Hamiltonian 
of the system which, however, is non-Hermitian since co- 
herent potentials are imaginary quantities. 

In order to solve Eq. (||) for the Green's functions, 
we will follow the procedure described below. First, 
let's solve the Shrodinger equation with the Hamiltonian 




(i = l,2). 



We choose these independent solutions so that two func- 
tions ipi{z) and ^pi{z) would describe two waves corre- 
sponding to the s-like electrons, and functions ip2{z) and 
(^52(2) would correspond to the d-like electrons. Namely, 
in a layer Fi (z < a) the solutions (pi{z) have the form 



<Pi(^) 



V2(^) 



exp[— ifcj'^z] 







z < a, 



z < a. 



(8) 



exp[— i/c^'^z] / 
and in a layer F3 (z > 6) the solutions ^i{z) are 



exp [jA:3*z] 




z > b, 



(9) 



f3=s,d 



Hf{z) #(z) = 0. 



(7) 



The solutions of this equation can be easily found since 
the potentials U^{z) have a step-like form. Let us put 
e = ep + where e^^ is the Fermi energy, and introduce 
the following notations: 



"Id 



are Fermi momenta in Fi and F3 ferromagnetic layers for 
electrons with the spin /i, and 



'i2s 







exp [i/cg'^z 



z > b. 



+0 (i 



1,2; a 



Since e = Ef + iO, then Imfcf" 
s,d). Thus these solutions satisfy the following boundary 
conditions: 

V'i(z)-^0 if z -> -1-00 (i = l,2), 
(pi{z) — > if z ^ —00 [i = 1,2). 

Starting from expressions (|^) and the solutions (piiz), 
ipiiz) can be easily extended in two other layers. Let us 
introduce the matrices 

$fz) = ( ^'("^^ ^ 
I ^f(z) ^^(Z) i ' 



Let also fcf " = ^ — and /cg" = \/ {^saY ~ '^^ 

be the components of electron momentum with spin fi 



The Wronskian of the system ^ is 
A = $^(z)L>^*(z) 



(10) 



where the matrix operator Dz is defined similar to Eq. 

, and the subscript T denotes the transposition opera- 
tion. It is known from the theory of differential equations 
that the matrix A is a constant since it satisfies the equa- 
dA{z) 



tion 



d 



0. Taking into account that e — ep + iO 



and Im E"^ < 0, the solution of Eq. (|^) for the retarded 
Green's function can be written in the matrix form as 



G«(z,z') = $(z)[A^]-il'^(z'), if z<z' 
G«(z,z') = *(2)A-i$^(z'), if z>z'. 



(11) 



To find the advanced Green's function, one has to put 
e = Ei? — iO in Eq. (||) and assume that Im E"'' > 0. 
Then we obtain 



(12) 



In the expressions (|11|) and ( |12[ ) the Green's functions 
depend on the in-plane momentum k and on the spin 
fi of the electron because solutions V'i(^) ^^id ipi{z) also 
depend on fi and n. 

Next, one has to find the two-point conductivity (H) 
of the system using the Green's functions ( pi] ) and (|12[). 
For that, we introduce the current matrices and 
(/i =1, 1) constructed with the use of solutions i^ii^) ^-nd 
<y9i(z), respectively: 



i^Hz)Dz^{z), 



(13) 



The total conductance of the system (T^p(z,z') may be 
presented in the usual form as a sum of the "bubble" 
part and the vertex corrections (see Fig. 2)l2J: 

a^p{z, z') = J^pfjO (z, z') + oil (^' + ^Mp(^' (14) 

The "bubble" contribution to the conductance de- 
scribes direct electron tunneling from the electrode to 
another one through the barrier when electron momen- 
tum parallel to the F/0 interface k = fcy is conserved 
(the specular transmission) . The scattering on the inter- 
faces affects direct tunneling so that the effective height 
of the potential barrier increases because the electron also 
has to pass through i5-like potentials at the points z = a 
and z = h formed by the self-energies S^(a) and 2^(6). 
The vertex corrections to the conductance describe the 
tunneling assisted by interfacial roughness — that is the 
processes when the electron with momentum k comes to 
the impurity center on the interface, undergoes scatter- 
ing (k — > k') in another channel with k' ^ k, and then 
goes away to the electrode or to the barrier. 

It also can be shown that the contribution to the 



tunneling conductance is negligibly small (~ e 



-49^(fc-a) 



from the diagram containing both vertex parts V a and Ffc 



a^p(z,z') = 




A^ P, A„ 



P, 



FIG. 2: The diagrammatic representation of the total two- 
point nonlocal conductivity ap,p{z,z'') as a smn of "bubble" 
and "vertex" parts. Here the full lines correspond to the 
Green's functions G^(z, z') and G^{z, z'), and wavy lines de- 
note the asymmetric gradient operator Vz of velocity at the 
points z and z' . The shaded square designates the vertex 
part r'^'' at the interface; h.c. denotes the complex conjugate 
terms. 



compared with other contributions that are of the order 

gram with vertices Fq and F^ 
corresponds to the interference of waves scattered from 
both interfaces. Therefore, we can neglect this interfer- 
ence term. 

Substituting obtained expressions (pT|), ([T^) for the 
Green's functions and using the definition ( |13| ) for the 
current matrices, we come to the following results. The 
analytical expression for the "bubble" conductance is 
given by 

(15) 

here it is assumed that z < z', A^ is the Wronskian (pT|), 
A denotes the junction area. 

For z < a, z' > b, the vertex corrections from the left 
and from the right interfaces can be written as 



'^I'piz.z') 



27rnA2 



-^ESp{a^«''^(^,^) (16) 

^ a;^?(&,^')} 



X F 



( Pi 71 

/32 72 



where summation is also performed over repeating indices 



Pi 71 



and Ft' 



are 



A,7. = s,d. HereF^^'^ 
the vertex parts on the interfaces a and b. A^„ are (2 x 2)- 



.PP ( Pi 71 
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matrices, their components are defined by 



(a)ImE«rn«)}- (19) 



Here the summation is also performed over repeating 
indices % and a^. The fulfillment of (^9|) provides the 
necessary condition of the nondivergence of the current 
through the tunnel junction: 



where [pfiaf^''' 
with the components 



and [p^(a; 



1 7172 



are density matrices 



K(«)] 
K(«)] 



0102 

ik 

7172 
ik 



s,d), 
s,d). 



(18) 



The expressions similar to expressions ( |l7|) may be also 
written for the matrices A^„(z, 6), Ap„{6, z') and for other 
position of the points z, z' with respect to a and b. The 
vertex parts and T'^'' are calculated in the "ladder" 
approximatiorO. The derivation of the equation for F^'' 
is presented below in the subsection IID. 

We have to note, that coherent potentials S^(a), S^(6) 
calculated in the framework of the CPA and vertices F^'', 
F^'' calculated in the " ladder" approximation satisfy the 
so-called Ward identity which in our case, for example, 
for the interface z — a can be written as follows (for the 
details, see Appendix A): 



(a) 



p=Ta 



-An 



Im G^^^Ha) 



(20) 



P=T,i 



dz' 



According to Eq. (|2C|) , the total conductance of the sys- 
tem is a constant value. In view of that, we will derive 
the exact expression for CT^p(z,z') evaluating the con- 
ductance at points z — a — and z' = & -I- 0, i.e. at 
the left and at the right sides from the interface. The 
Wronskian matrix and the current matrices j'^{z) 
and j^{z') are expressed in terms of matrices $(z) and 
^'(z) — these matrices, as follows from Eq. (pi]), deter- 
mine the Green's function. The straightforward evalu- 
ation of the conductance according to formulae (|l^-[r^) 
leads to the result that a^p{a,b) is expressed in terms of 
the retarded Green's functions at the points of interfaces 
Gp{a) = G^iz ^ z' = a) and G.(6) ^ Gj}{z ^ z' ^ b) 
constructed according to Eq. (11). The explicit form of 
Gp{a) in the (k, z) representation is given by the expres- 
sion: 



2™, 



■(* 



cotan (pi^) 



fi d 



2mA 



-s^'(«) 

cotan {kfzQ + (fi"^) 



(21) 



where 



us fci ITT-^ 

t&nip^ = 



tanp'^ = -1 

92 "^rf 



The similar expression can be written for the Green's 
function Gp(6) on the right interface. For that, one has 
to make the substitutions fc^ fc^", ^'j^>^{a) S2^(6), 
93^ " — > (^3 " in the expression for G^ (a) . Here ttIq, (a — 
s, d) and jti^j.^^ are the effective electron (hole) masses in 
the ferromagnetic and insulating layers; fc^*" and ^2 (* = 
1,3, a — s,d) are the functions on n introduced above 
in the text after Eq. (^. Let us define the "transport" 
density of states as follows: 



A,{b) 



Glia)j^G^{a), 
Giibfj^G.ib), 



(22) 



where 



■J h 







\ 







^ d 



_3_ 







^ 

rrid / \ md 

Expression ( p^ ) for the "bubble" conductance then 
reads 

„2 



2'KnA 



where 



e«2^o 
A" 



xAlia)X-'[^)x^'A,ib)j, (23) 
. _ sin(/cf'^zo-f (^5"*) 



11 d 

sm ipi 
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£92 20 





\ — 



a d 



( sL 



\ 



-92™ 



92 

m2 



-«2' 



and w — b ~ a ~ 2zq is the "width" of d-barrier. 

For the vertex correction from the left interface (z = a) 
we obtain 



Ap(a,a) 
Ap(a,5) 



xrr(^: (24) 

Gt(a)A-i(^)A,-Mp(6) 
X A,-i(^)A-iGp(a). 



In the similar way, the vertex correction from the right 
interface (z — b) reads 

A,ia,b) = Gt(6)A-i(£)A-M^(a) 

X A-^(£)a,-1G,(6), 
Ap(6,6) = -Ap(&). 



Expressions (^3|-g5|) determine the total conductance of 
the system. 



C. The CPA equations 

The scattering of the electron by random potentials 
within the interface (terms -ffspin-cons. and iygpin-flip hi 
the Hamiltonian) is taken into account by introducing the 
self-energy operators EJ^''(a) and into the effective 

single-particle Hamiltionian (|^). To calculate the self- 
energies, we apply the coherent-potential approximation 
(CPA)HM Let's denote kets \j,fi,pZ) = \pZ) 
corresponding to the Wannier states of the electron on 
the interface a (a — 1,2) at the given site in the 
(x, ?/)-plane, where 7 refers to the state s or d and /i is 
the spin. The symbol (E) denotes the direct product. The 
problem of finding of the single-particle Green's function 
G^p(r,r') of the many-body Hamiltonian is reduced 



to the related single electron problem of the propagation 
of the electron in a random interfacial potential 



n:Q = l,2 



where v" — Vn^^" + Vn''" and the summation goes over 
the interface number a and over the sites n. Here 

i'i"""^ E ln{\s,^^){d,p\ + \d,^i){s,f,\} (26) 

is the random potential of s — d hybridization associated 
with the site n and 

^" = Jn {k, T)^-(p")(s, 1 1 + l-s, i)s+{pZ){s, T 1} 

(27) 

is the exchange-like interaction with the surface magnons. 
The random quantities 7" and J" used here were intro- 
duced in Sec. HA. 

Now one can formulate the CPA by the ordinary way 
and the only difference with respect to the usual situation 
of the bulk scattering is that the initial Green's functions 
have to be calculated for the trilayer system. We assume 
that the averaged propagator of the system G|^''(r, r') 
differs from the initial Green's function, corresponding 
to the Hamiltonian Hq [see Eq. (|l|)] , by the self-energy 
correction in Eqs. (|^), (^. This means that the system 
behaves as if coherent potentials E^'' had been assigned 
to each site of the interface a and b. After the introducing 
of the effective medium , each site becomes a source 
of the random potential = v" — S". The single-site 
t-matrix associated with potential is given by 



~,(sf)c 



1 



1 - {< - E"}G(z„) 
where = a or b, 



K - S"}, (28) 



= E E IAM>S^"(^a)(7,Ml, (29) 
= E E l/3,M)G^^(^a)(7,Ml- 



Here 



(30) 



is the averaged Green's function at the interface a which 
is expressed via in accordance with Eq. (21). The 
upper limit Kmax is a cut-off of the in-plane momen- 
tum which originates from the finite size of the Bril- 
louin zone. For that, one has to substitute the Brillouin 
zone's projection onto (A:^., fcy)-plane by a circle of the 
same square with a radius Kmax- For the bcc lattice. 
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i^max = 2Y7r/\/2/ao7 where ag is a lattice constant. The 
single-site i-matrix ( [28| ) is obviously diiTerent for differ- 
ent sites. At the same time it is supposed that scattering 
by the random potential is taken into account in the av- 
eraged propagator G{za) by the self-energy operator S". 
Therefore, we require that the ensemble average of the 
single-site t-matrix vanishes, i.e., 



Alb ■ 



Bib 



0. 



(31) 



Here and t'g are the values of the single-site i-matrix 
associated with a given site n which is occupied by the A- 
type ion (ferromagnet's ion) or by the B-type ion (impu- 
rity), respectively. Brackets denote the averaging 
over magnon degrees of freedom. The equation ( |3l| ) is the 
well-known selfrconsistent coherent potential approxima- 
tion (CPA)E3'E3 that implicitly determines the self-energy 
operator S". The CPA equations are formulated for 
the particular case of electron scattering within the F/0 
interface. The straightforward calculations show that the 
single-site t-matrix t" can be represented in the form: 



(32) 



with respect to spin-up and spin-down subspaces. Here 
71+ = S+S^, n_ = S-Sj^, and the blocks f"^, iZ^''^^ are 
(2 X 2)-matrices, functionally depending on n+ and n_, 
with the components designated by indices of bands s 
and d. 

To satisfy the condition (^l|), one has to consider 
only the spin-conserving part of this equation as long 
as (t"^S'+)b = since the expression to be averaged con- 
tains unequal number of the creation and annihilation 

operators of the magnons. In order to calculate (in^''^'')fc 
we adopted the further approximation and assumed that 

{tZHn+))b = iZHn), 
rJ{n-))b = iZHn), 

here n = n{T) is the average number of magnons at the 
given temperature. In other words, we substituted the 
operators h± by its averaged values. The function n{T) 
is given by the familiar expression 



niT) 



d2q 



1 



1 



(33) 



where Lu(q) is the spectrum of surface magnons. The 
approximation we made to take into account the spin- 
flip processes is the static approximation for magnons. 
The inelasticity of the electron-magnon processes may 
be taken into account in the energy conservation rule. 
However, the characteristic magnon energy htoo is much 
less than the Fermi energy ep, and in the first approx- 
imation one can neglect this energy Hloq — it becomes 
important in the case of finite voltages when the process 
of emitting of a magnon by a hot electron influences on 




T -\- T 



"2 




«2 yi 



FIG. 3: The diagrammatic representation of the equation (see 
the text) corresponding to the calculation of the vertex part 
r'''' in the "ladder" approximation. 



the form of / — dependence. We, therefore, restricted 
our calculation to the case of small voltage bias with the 
voltage less than huiQ. Within this approximation, for 
the system of CPA equations we get 



e^{n)x + tl^{n)y = 0, 
t'XHn)x + t%^{n)y = 0. 



(34) 



As far as matrices tZ^''^^ are (2 x 2)-blocks, thus for the 
general case Eq. (|^) represents the system of two matrix 
equations for 8 unknown quantities E^'''(zq,). The system 
( p4| ) is one of the possible forms of the CPA equations. 
But actually we used another representation of the CPA 
which was more convenient for the numerical implemen- 
tation. For that we ejiploited the augmcntcd-spacc for- 
malism by MookerjeeE3, and the details are presented in 
Appendix B. 



D. The vertex corrections 

To find the vertex corrections we used the "ladder" 
approximation combined with the CPA, that was orig- 
inally proposed by Velickya. The diagrammatic rep- 
resentation of the "ladder" approximation is given in 
Fig. 3. Since the scattering potentials on the differ- 
ent interfaces (at the points a and b) do not correlate 

with each other, the vertex parts T^f,-, | '^^ 

be found independently for each interface and we will 
omit interface suffix in the subsequent expressions. Let 

01 Ji \ _ /(-J.Q^*/3l7l /j-a-i 72/32 

the product of two random i-matrices over configurations 
and magnon distribution, where as before /i, p are spin 
indices and (3, 7 are orbital indices. On the diagram in 
Fig. 3, it is presented by the dashed line. The analytical 



can 



be the average of 



expression for the correlator T^'' 



can be found 



in accordance with the adopted approximate scheme of 
averaging over the magnon degrees of freedom. Taking 
into account Eq. (|32|) we obtain: 



j-TT ( ft 71 
P2 72 



*| , xft7l I , .72/32 
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72 02 
72/32 



(35) 



ynt 



\^'i7i , 



01 71 

02 72 



a;nt^ (n) 



/3l7l,+ / \72/3: 



72/32 



t ( W -4- 

where t2(^^^{n) and t^^^^jn) are the components of the 
single-site t-matrix, Eq. (p2[). We also define the operator 



- G:'3i7i(^^)Gfe72(^j 



(36) 



denoting the propagator of the pair of electrons in the 
"ladder" equation at the interface Zq.. Its definition fol- 
lows from the fact that in the diagram representation 
of the CPA the multiple scattering on the given site is 
assumed to be incorporated into the single-site i-matrix 
i", corresponding to the single vertex of any diagram. 
Due to that, the subsequent sites in the "ladder" di- 
agrammatic equation (Fig. 3) must not reproduce each 
other. Therefore, the necessary correction is subtracted 
in Eq. (36). The summation over k goes up to K^ax 
similar to Eq. (30). After that definitions the analytical 
equation for the vertex part reads as 



YfJ-P I "1 /3i 

a2 02 



ai 01 
a2 02 

71 <5l 



- Tfj,\ 

pAp 



ai 71 



<5l 01 
S2 02 



(37) 



Eq. (|37| ) is the ordinary system of linear equations. It can 
be proved that the vertex parts found with the use of the 
"ladder" approximation and self-energies found with the 
use of the CPA satisfy the Ward identity (|l^) , and hence 
the total tunneling current is a constant value j{z) — jo, 
i.e. it does not depend on z. The proof of the Ward 
identity is presented in Appendix A. 



III. RESULTS AND DISCUSSION 

We considered the case of Fe/Al203/Fe tunnel junc- 
tion. The following parameters were chosen to describe 
the system. According to estimations of StearnscJ, for 
the itinerant s~like electrons in Fe we set /cf ^ = 1.09 A~^, 
= 0.42 A~^, and = 1.0 me (here rUe is bare elec- 
tron mass). For the more localized d-electrons (holes), 
we put rrici ~ — lO.Owe. The d density of states r-at 
the Fermi energy is larger for the minority spin bandO, 
and we put = 0.45 A~\ 



k^^ = 1.35 . The values of the Fermi momenta define 
the positions of band bottoms V/^, V^^ {a = s,d, fj, =t, i) 
with respect to Ef- For d-electrons the Fermi momenta 
were chosen such a way that the interface density of 
states Ps(j;)(e_F) = ~1/^I™G'^j-^j(zq,) (where Za = a,b 
are the positions of interfaces) comply with the follow- 
ing ratio pi^^Hep) : pliep) : piisp) ~ 0.1 : 1 : 10, 
which is the typical situation for the case of 3d transition 
metals (see, r&r a example, the calculations of Tsymbal 
and Pettifoio). We also put the width of the barrier 
w = 20 A, and zq = ao/2 where ao = 2.87 A is a lat- 
tice constant for bcc Fe (zq is a parameter describing the 
thickness of the interface). 

The main features of the band structure of a-Al203 
crystals were briefly presented in Sec. IIA. In a view of 
that discussion, the following parameters of the model 
were taken to describe the amorphous AI2O3 barrier: the 
effective masses of electrons and holes are = 0.4 me, 
= —10.0 me, the positions of the conduction band 
bottom (Us) and the top of a valence band (Ud) are Us = 
— Ud = 3.0 eV, i.e. the width of the band gap is ~ 6.0 eV, 
and is assumed to be a zero of energy. 

In order to illustrate the general formalism presented 
in the previous sections, let us first consider the case 
of only s — d impurity scattering (i.e. T = and there 
are no spin-flip processes) when concentration of impu- 
rity ions on the interface is (1 — x) = 0.5. In this case 
all formulae have a simple analytical form. As it was 
mentioned previously in Sec. IIA, in the two-band tight- 
binding model description of a binary alloy the random 
variables are diagonal matrix elements of the Hamilto- 
nian £^(cc), £^(x) and jn{x) depending on whether the 
site n is occupied by an A or B ion. To be more precise, 
one must consider the concentration dependence of these 
matrix elements. It reflects the fact that their values are 
modified by the existence of a charge transfer in the al- 
loy. Since in the framework of our model we are focusing 
on the hybridization effects (see Sec. IIA for justification 
of this approach), we have 7^(2;) ^ Jq, Jb{x) 7^ 7(f 
where and Jq are hybridizations of the pure com- 



ponents. Let us assume for simplicity that in Eq. (B2) 
70 = 7 = xjA + UlB = 0. Then for x — 0.5 we get 
jA — —7s, and thus 7^ = x{l~x){'yA — ^B)'^ is the single 
parameter characterizing the amplitude of s-d scattering 
on the interface. The nonzero value of 7 will lead to a 
slight modification of the s and d-like eigenstates describ- 
ing s and d-like electrons in the vicinity of the interface, 
and thus will not affect the qualitative results presented 
below. We will omit spin suffixes since we consider the 
spin-conserving scattering now. At x = 0.5 only diago- 
nal elements of the self-energy matrix T,"^ in the sd-space 
have non-zero values. Then the Green's functions, e.g. at 
the point a, become [see Eq. (^) ] 



F I 

therefore, k^ > k. 



FT 
d ' 



k^ 
2m. 



{i — cotancpj) 
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(38) 



L 2mrf 

-1 

dd 



[i - cotan {kfzo + Lpf)] 
kcIk 



and CPA equations are written as follows 
^'G^' ^dd ^ l^G^ 



1 + Y^ddQdd ' 



1 + S^'^G" 



(39) 



which must be solved self-consistently by means of con- 
verging iterative procedure. As far as Green functions 

are diagonal, one can denote = F ^ " |^ ^ and can 

introduce the similar notations for quantities T and -D, 
defined in section IID. Then one gets 



T, 



7 



7 



\l + Y.ddQdd\'' |1 

Tss = Tdd = 0, 



1 - T^.D^Wdd' 

sd 



dd 



(40) 



1 - T^^D'^'Ddd' 



1 - T'^^D'^D'^'^' 



Since the mass of holes 



in the insulator is much 



larger than the electron mass to^, the exponential factor 
g-2i?2(&-a) formulae (^ - ^) for tunneling conduc- 
tance is negligibly small as compared with one for s-like 
electrons. Therefore, one may neglect the contribution 
from d-holes to the tunneling current. Then for the tun- 
neling conductance we obtain 



where 



e 

2Tih 



'^^ ^ (5)e-29>^ 



(41) 

is the conductance corresponding to the "bubble" dia- 
gram (here w = b — a), 





2'Kh 



-2q2'W 



(42) 
27r 



27r 



is the "vertex" contribution to the conductance on the 
left interface, and 



2nh 



ndn 



2 

2/'_92_ ' 

m5 



(43) 

n'dn' 
27r 



27r?i 



+ 7^ / ^r.W^E r(.^^ 



27r 



|2/' 92 ^ ASS /„\„-2qt,w 1^' 



2tt 



is the "vertex" contribution to the conductance on the 
right interface. Here are the transport densities of 
states which for the present case on the left interface are 
given by the expressions 



\kl[i — cotan c/?|) /2ms — S^''('j)| 

\kf[i - cotan(fcfzo + (pf)]/2md - T.df{a) 



(44) 
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and by the analogous expressions in the case of right 
interface. 

We remind, that physically the "bubble" term in the 
total conductance is a contribution to the current due to 
direct tunneling when the electron momentum k = fcy in 
the plane of the interface is conserved. Scattering pro- 
cesses renormalize the " bubble" term with respect to the 
case of ballistic transport such that the effective " height" 
of the potential barrier seen by electron increases due to 
self-energy corrections E^* arising on the interfaces. The 
vertex contributions to the conductance describe the tun- 
neling assisted by interfacial scattering — that is the pro- 
cesses of tunneling with scattering on the F/0 interface 
when the in-plane momentum is not conserved, i.e. k' ^ k 
for the scattered electron. Thus, both momenta k and n' 
in Eq. (|6|) determine the vertex corrections. 

Note also, that both vertex corrections ^-nd ar^ 
( ^ ) , (IH) consist of two terms with vertices Tgd and . 
The terms with vertex part Tss are contribution to the 
tunneling conductance from s-like electrons only. The 
scattering in the s channel described by (^) is caused 
initially by s-d scattering. The terms with Tsd describe 
either the process of diffuse scattering of d-like electron 
to the s state in the F/0 interface and then tunneling of 
s-like electron in the barrier, or the process of tunneling 
of s-electron and then its scattering to d state in the 0/F 
interface and leaving into the electrode. 

We have calculated the tunneling conductance and the 
TMR ratio defined as TMR = (cp — crAp)/o'AP, where 
CTp and (Jap are the total conductances for the paral- 
lel (P) and antiparallel (AP) alignment of magnetic mo- 
ments in the F-layers. The CPA equations defining the 
self-energies were solved numerically. The validity of the 
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FIG. 4: "Bubble" and vertex contributions to the tunneling 
conductance (in the units /2-Kh per unit square lA^) for the 
parallel (P) alignment of magnetic moments in the ferromag- 
netic layers as a function of the scattering parameter 7 on the 
interface in the absence of spin-flip processes. The param- 
eters of the model are: s-like electrons — " = 1.09A~^, 
fcf ^ = 0.42A-\ = 1.0, w}l = 0.4; d-like holes — 
fcf = 0.5A-\ fcf = 1.4A-\ ■md = ml = -10.0, the height 
of the potential barrier Ua ~ —Ud ~ 3.0 eV, the width of the 
barrier w — 20 A, the concentration of Fe ions on the interface 
X = 0.5. 



Ward identity, Eq. (|l|), was checked after the vertex 
parts r had been computed at every step for a given 
value of the parameter 7. 

The results are presented in Figs. 4-7. In Fig. 4 the 
"bubble" (|l|) and the "vertex" (||), @ contributions 
to the conductance for spin-up and spin-down channels 
are shown for the parallel (P) alignment of magnetic mo- 
ments as a function of the scattering parameter 7 on 
the interface. (We remind that up to now only the spin- 
conserving scattering is considered.) As it was mentioned 
previously, scattering suppresses the "bubble" conduc- 
tance, and hence cr°| and a^^ are decreasing functions 
on 7. The contribution from the "bubble" part is larger 
for the majority spin (t) channel since fej^** > k^'^. On 
the contrary, the contribution from the vertex corrections 
dominates for the minority spin (|) channel — that be- 
havior can be understood as follows. 

First, at small 7^ the imaginary part of the self-energy 
S'"*, which describes the scattering of s-electrons, be- 
haves as ImE^jj^j ^ 7r7^/9^^j(ei?) [see Eq. (|9|)], where 
p^^^{eF) are spin-up (down) d density of states on the 
interface. As far as pfisp) is of the order of magnitude 
greater than p'^{eF), the scattering for s-like itinerant 
electrons from minority spin band (J,) is more effective 
than for the majority spin (t) electrons. The more effec- 
tive scattering of spin-down s-electons leads to the pre- 
dominance of the minority spin contribution with vertex 
part r|j, over the majority spin contribution with vertex 



FIG. 5: "Bubble" and vertex contributions to the tunneling 
conductance (in the units /2-Kh per unit square lA^) for 
the antiparallel (AP) alignment of magnetic moments as a 
function of the scattering parameter 7 on the interface in the 
absence of spin-flip processes. Parameters of the model are 
the same as in Fig. 4. 



Tl^, since at small 7, tI\^> ~ T^^fu) ^ ' 

therefore F^j, ^ Fj^ Second, the transport density of 
states A'^'^{a) [Eq. (p4|)] which determines the partial con- 
ductance due to mixing of s and d channels in the in- 
terface (terms with T gj) is also larger for minority spin 
channel. Thus, the combination of these factors results in 
domination of the vertex contribution to the conductance 
from the minority spin channel. The results of calcula- 
tions also showed that the ss contribution to the vertex 
corrections (terms with vertices T ss) was much less 
than the sd contribution a^^ (terms with Vsd)- 

The " bubble" and vertex contributions to the conduc- 
tance (which are the same for spin-up and spin-down 
channels) for the antiparallel (AP) alignment of moments 
as a function of 7 are presented in Fig. 5. The total 
conductances for P and AP configurations are shown in 
Fig. 6. Finally, in Fig. 7 the TMR ratio (ctp - aAp)/(yAP 
as a function of scattering parameter 7 is shown. When 
the amplitude of scattering is negligibly small (7 = 0), 
we have positive value of the [XMR ~ 24% — that is the 
result of Slonczewski's theoryll3 under the chosen param- 
eters for s-like electrons. With increasing of 7 the TMR 
amplitude is monotonically decreasing and can become 
even negative if 7 > 7c ~ 0.16 eV. 

To understand qualitatively the obtained results, con- 
ductance and magnetoresistance can be rewritten in the 
general approximate form: 

dp ^ {p^ ) + [p^ j + P3 1 ^^Pi + P3 ^ sdPl ' 
(JAP W ^Pl Pi + Pz^ sdPl + P?. KdPl ^ 

Act — ap — gap (45) 

« (pf - Pi'? + (pf - pf )(pf f L - pff L), 
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FIG. 6: The tunneling conductances of the individual spin 
channels (in the units /2-Kh per unit square lA^) for the 
parallel (P) and antiparallel (AP) alignment of magnetic mo- 
ments as the function of the scattering parameter 7 on the 
interface in the absence of spin-flip processes. The parame- 
ters are the same as in Fig. 4. 



where we took into account that <C crf'^, and for high 
and thick enough barrier the main contribution to the 
tunneling conductance is due to electrons with momen- 
tum almost perpendicular to the barrier (factors e"^*^™). 

Consequently, p^^'''^'' are the quasi-one-dimensional s den- 
sity of states near the interface, and ^3^''''^'' are the cor- 
responding three-dimensional d density of states. The F 
are renormalized vertex corrections. In expression ( |45| ) 
for Act, the first term due to direct tunneling is always 



positive. Concerning the second term, > 



> k 



si 



but 



> Pa 



si 

p^' smce 
and, therefore, this term is 



negative. This contribution decreases the magnetoresis- 
tance when scattering parameter 7 becomes larger. At 
7 > 7c the sd contribution overlaps the contribution to 
Act from the ss channel, and we have TMR < 0. Thus, 
the inverse TMR ratio arisen in our model is caused by 
extremely strong scattering of negatively polarized d-like 
electrons (which give the indirect contribution to the tun- 
neling current) to the s-band on the interface. 

The parameter 7 determining the amplitude of s-d 
scattering on the interface is defined as 7^ — x{l—x) (7^ — 
7b)^- One can regard the difference between hybridiza- 
tions, 7^1 — 7_B, as approximately a constant value for 
given constituents. However, (1 — a;) is a concentration 
of impurity centers, and thus 7 is a measure of the imper- 
fection of the F/0 interface. The proposed model, there- 
fore, explains qualitatively the strong intesface sensitivity 
of the tunneling magnetoresistance effectQuB. According 
to Heiiie's discussion of the hybridization in transition 
metalsO, the hybridization constants are from about 1.0 
to 3.0 eV for different elements. One can assume that the 
difference (7^ — 7b) is of the order of magnitude smaller 
and, therefore, the critical value of scattering parameter 
7c ~ 0.16 eV that we obtained is more or less realistic. 



FIG. 7: The tunneling magnetoresistance (TMR) as a func- 
tion of the scattering parameter 7 in the case of only elastic 
scattering on the interfaces. The parameters are the same as 
in Fig. 4. 



The strong reduction of the TMR due to nonideal 
structure of the metal/insulator interfaces is a well- 
known observation. The oxidation of a thin Al layer 
leads to the undesirable oxidation of few metal mono- 
layers close to the F/0 interface and thus to the. for- 
mation of F-0 oxides ^fep04ll3, CoO and C03O4B) — 
that reduces the TMRilli On the other hand, if the 
too thick Al layer is not oxidized completely, the con- 
tamiiiatipn of the interface by Al ions also reduces the 
TMrMeI. The dependence of the TMR vs. the thickness 
of the Al overlayer has, therefore, a maximum and the 
best TMR values-achieved by Moodera's group lie in the 
range 10 16 M. 

The contamination of the F/0 interface by OH ions in 
the early experiments by MerserMey and Tedrow on tun- 
neling with the superconductoril led to the small mea- 
sured values of a spin polarization P for Ni and Gd. The 
contamination was due to oxidizing of the Al films in 
the laboratory air. The improved technique of samples 
preparation in a pure oxygen increased the vaJHCs of P 
for Ni and Gd, and for some rare-earth metalgSfl. 

In a recent work by LeClair et aZHj the strong suppres- 
sion of magnetoresistance was observed in C0/AI2O3/C0 
tunnel junctions with a very thin Cr interfacial layer. It 
was found that the TMR decayed exponentially on the 
Cr interlayer thickness with a length scale ~ lA (approx- 
imately 0.5 monolayers). With the addition of only 3 A 
Cr (« 1.5 monolayers) the reduced TMR was only 10^ 
of the initial value for a control junction. LeClair et a/.Ej 
presented some qualitative arguments that the suppres- 
sion of a spin polarization (and, hence, the reduction of 
TMR) was due to more strongly suppression of major- 
ity s-p density of states compared with minority spins 
caused by the resonant scattering of majority spin sp- 
electrons with the Cr d states. From the point of view of 
our model, we can explain the strong degradation of the 
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FIG. 8: The TMR for different temperatures as a function of 
the scattering parameter 7. 



TMR by the strong electron scattering within the inter- 
facial Co-Cr alloy that is formed under the preparation 
of extremely thin (~ 1 monolayer) Cr interlayer. 

We have also calculated the temperature dependence 
of the TMR taking into account spin-flip scattering in 
addition to s — d impurity scattering as it is described 
in details in sections IIC and IID. For that, the average 
magnon number n{T) ( |3^ ) as a function of the temper- 
ature was found in analogy with Debye's treatment of 
phonons using the-, similar approach that was proposed 
by S. Zhang et a£3. The magnon dispersion relation in 
Eq. (^3|) was replaced by simple isotropic parabolic spec- 
trum 



— En 



Curie temperature Tc and in the mean-field approxima- 
tion is given by — ?>kBTc/{S + 1). For the chosen 
model of the dispersion relation lUq, one has to over- 
come the divergence on the lower limit of the integral 
in Eq. (^3|). Therefore, one must introduce a lower wave- 
length cutoff Ejt^. Physically, it may represent a finite 
coherence length due to interfacial roughness. In our cal- 
culation we have takeii-|the same parameters that were 
used by S. Zhang et alx^ for the analysis of the zero-bias 
anomaly: S = 3/2, keT^ = 110 meV, and ^ ^ meV. 
Then, for the temperature range within the room value, 
we have 



n{T) 



-EJksT 



where /cmax is the equivalent radius of the two-dimensi- 
onal Brillouin zone (see Eq. 30), and Em is related to 



We have put x — 50%, Ja = 2.0 eV (for Fe atoms) 
and Jb = eV, i.e. it is supposed that spin-flip process 
is possible, if an electron scatters on the Fe ion. The 
TMR vs. the scattering parameter 7 at different temper- 
atures T = 4.2, 77, 210 and 300 K are plotted in Fig. 8. 
In Fig. 9 the temperature dependence of the resistance 
for the P and AP configuration of magnetizations is pre- 
sented for the same parameters, but 7 = 0. The results 
show that the TMR ratio decreases with increasing of the 
temperature. Moreover, the resistances Rp and Rap of 
the structure for both configurations, P and AP, are also 
decreasing when temperature increases — that is in the 
qualitative agreement with experimental data (e.g., see 
Fig. 4b in Ref. [38]). The physical mechanism of this ef- 
fect is related with the excitation of spin-flip processes in 
the system. Due to these processes, the new channels of 
electron scattering appear which are frozen at zero tem- 
perature. As the result, the conductance of the system 
increases for both P and AP configurations and, there- 
fore, the resistance drops. The spin-flip processes mix 
the spin-up and spin-down channels. Therefore, the rel- 
ative difference of the resistances decreases at different 
configurations and the TMR decreases with increasing of 
temperature. 



0.44 




0.20 



50 100 150 200 250 300 

T(K) 



FIG. 9: The temperature dependence of the resistance for P 
and AP configurations of the magnetic moments at 7 = and 
X = 50%. 



IV. CONCLUSIONS 

In conclusion, we would like to outline the main re- 
sults obtained in the present work. Based on the analy- 
sis of the band structure of 3d ferromagnetic metals and 
AI2O3 crystals and on the results of ab-initio calculations 
of the magnetoresistance for epitaxial tunnel junctions, 
we built on an adaptation of the simplified two-band s-d 
model to treat the diffuse electron transport in the non- 
ideal F/O/F magnetic tunnel junctions. We had mod- 
eled the rough F/0 interface by the random binary alloy 
that is formed from the ions of the ferromagnet and im- 
purities (e.g. the Al or O ions), and assumed that the 
main mechanism of electron scattering on the interface 
which substantially affects tunneling is the s-d scatter- 
ing. We used the Kubo formalism to calculate the tun- 
neling conductance and found the vertex corrections to 
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conductivity with the use of the " ladder" approximation 
combined with the CPA. The obtained results show that 
in the case of strong electron scattering within the in- 
terfacial alloy the vertex corrections give the essential 
contribution to the tunneling conductance. We proved 
that adopted approximations lead to the physically cor- 
rect results, namely, the non-local conductivity tensor is 
a constant function and, therefore, a tunneling current is 
conserved. We showed that interfacial inter-band scat- 
tering substantially reduces a value of the TMR, which 
can become even negative in the case of extremely strong 
scattering. The reason of the suppression of the mag- 
netoresistance is the indirect contribution of negatively 
polarized d-like electrons to the tunneling current due to 
strong scattering to the s band on the interface. It is 
also shown that spin-flip electron scattering on the sur- 
face magnons within the interface leads to a further de- 
crease of the TMR at finite temperature. Thus, in the 
framework of the proposed model, we are able to explain 
qualitatively the strong interface sensitivity of the tun- 
neling properties that is observed in the experiments. 
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APPENDIX A: WARD IDENTITY 

In this Appendix we briefly describe how to obtain 
Eq. (|9|) (Ward identity), and using a simple example we 
will show how this identity can be proved. 

Starting from Eq. (po|), one first has to compute the 
derivation of the "bubble" conductivity. Assuming, that 
the self-energy S^'^ is symmetric with respect to the re- 
arrangement of band indices, i.e. S^'^ = S^": a-nd using 
the fact that functions ip°'{z') are the solutions of the 
Schrodinger equation (^, for the derivation of the cur- 
rent matrix we obtain: 



d 



E 

a.f3—s,d 



where 



Ef (z') = J:f{a)S{z' -a) + ^f{h)5{z' - 6), 

and the matrix j£j^(z')]"^ is defined by Eq. (^. Then, 
from Eq. (15), (0) we have 



2-KhA 



^ ^ [2ImEf (a) 



X 5(z'-a)A«f(z,a) (Al) 
+ 2ImEf (6)5(z'-fe)A^f(z,6)\ 

To proceed further, let us compute the derivation of the 
vertex correction to conductivity. The initial expression 
for the matrix A^^ (see diagram in Fig. 2) is written as 



2ima 



Taking into account that the Green's functions are the 
solution of the Eq. (^), we obtain 

±.Kl^^^{a,z') = -25{z'-a)lTnGl^:^{a) 

+ 2 J2 ImEf (a)G;^^^(a) 

a/3— s,d 

X GP^l^ia)Siz'-a). 

Substituting the obtained expression to formula ( p^ for 
the vertex correction, we get 



27r?iA2 



Y,2Siz'-a)A''^tiz,, 



/3i 71 
^2 72 



ImG;^?(a) (A2) 



{a)G;:r^ia)G';i7'ia) 



The summation here is also performed over indices a^, 
f3i, and 7^. A similar expression can be written for the 
derivation of the vertex correction a^.'iiz, z') at interface 

I 1 A k 

z = b. From ( |A1| ) and (|A2[) one can obtain the final 
expression ( p9| ) for Ward identity. 

Let's prove identity for the simple case of only s — 
d scattering, when there are no spin-flip processes, and 
X = 0.5. This situation was considered in Sec. Ill, and 
there were introduced the notations Fq^, T^p, and D"^ 
{a(3 — s, d) for the components of the (2 x 2)~-matricies 
r, T, and D. Within these notations, the identity (p^), 
which has to be proved, can be written as 



l3=s,d 



IraG^^ --y\G^/flmE^^ 



Here we omitted spin suffixes, all values associated either 
with the interface a or 6, G^^ and Tap are defined by 
Exps. (pi Mj. From (M it follows that 



D 



\G° 



T-^ - D. 



K a()=s,d 



Using these expressions. Ward identity can be written in 
the form: 

ImG""- {[r-i]„/3 + 5„/3|G""p}lmE'5'' = 0. 

I3=s,d 
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Using (pol), we then can find that 
1 



Den 



'en 



Im^dd ^ _^(\Gss^2^^(jdd 

Den 



where 



Den = \G"'\'^\G'^'^\'^ - j'^'ll + J^^'G" 



The same expessions for ImE*'' and ImS'''^ follow di- 
rectly from the CPA equations ( p9| ) — that finishes the 
proof. 



APPENDIX B: THE DERIVATION OF THE CPA 
EQUATIONS 

In this Appendix we will derive the CPA equations for 
our particular case using the augmented-space formal- 
ism (ASF)Ea. As was described previously, we assume 
the F/0 interface to be a random binary alloy of the 
type AxBi-x, where A are ions of the ferromagnet and 
B are impurities. Following the ASF, we associate each 
random variable 7" and with the self-conjugate op- 
erators 7 and J, respectively, which are determined in 
the auxiliary 2-dimensional vector space $ such a way, 
that the spectrum of these operators coincides with the 
set of possible values of random variables. For the sake 
of clarity, hereafter, the sign "tilda" is ascribed to any 
operator acting on the auxiliary space. We also define 
the orthonormalized basis |s), where s = A or _B, which 
are eigenvectors of 7 and J, so that 



j\A)=-fA\A), 
^\B)=^b\B), 



J\A)^Ja\A), 
J\B)=Jb\B). 



(Bl) 



According to that definition, 7 and J commutate with 
each other. Let now 7(7", J") be a function or an oper- 
ator of random variables 7" and J". Then, the operator 
in the auxiliary space associated with function /, is de- 
fined as / = /(7, J) and according to (Bl), e.g. 



is a value of /, if the site n is occupied by an ion A. One 
can introduce another orthonormal basis in $ 

\Q) = ^\A)+^\B), \l)^^\A)~^\B), 

so that the operators 7 and J in this representation are 
written as 



70 7 
7 71 



J 




(B2) 



Here x and y are the concentration of A-type ions (fer- 
romagnet's ions) and B-type ions (impurities) on the in- 
terface, respectively, and 

70 = 2:7^ + ?;7_B, 71 2/7A + 2^75, 
7 = ^iy(7A - 7s), Jq = xJa+vJb, 
Ji = yJA + xJB, 5 = ^/xy{JA~ Jb)- 

Then one can prove that the averag e v alue of / is given 
by / = (0|/|0). Together with Eq. (|b|) this property is 
the way to evaluate the average of any given operator, 
depending on the random variables, and we can apply 
this method to average the i-matrix ( p8| ) . 

Following the general scheme of the ASF, outlined 
above, the random effective potential = — is as- 
sociated with the operator U —Y, acting in the augmented 
vector space $ (g) L, where L denotes the 4-dimensional 
space of orbital (s, d) and spin ("f , |) electron degrees of 
freedom. In accordance with (p6|), ( [27| ) and ( p32| ) this 
operator has the form 



U -Y. 



Uo-Y A \ 
A Ui~t j 

7o -I- Jo^' — S ^ + 5v 
+ 5v 71 + Jiv — S 



(B3) 



where S is defined by (p9|), operators 7(j) are defined 
similar to (Eg), and other operators are given by 



« - IT)5-(p")ai + U)^+(A")(TI, 

5 - 5\s){s\, J,^J,\s){s\ (i-0,1). 

Let us also introduce the nonrandom averaged propa- 
gator acting in the augmented space 




G = 



and associated with potential U 
scattering i-matrix 



t 



1 



GiU-Y. 



) the augmented 



(B4) 




Its projection onto the zero- level |0) of the augmented 
space ioo = (0|i|0) coincides with the average from the 
"physical" random t-matrix (]2^). The subsequent aver- 
aging over magnon degrees of freedom (ioo)b must vanish 
due to condition (|4|). 

To proceed further, let's introduce the electron propa- 
gator Gi 



Gi 



1 



G-i - (C/i - E) 



1 1 



which is associated with the "propagation" of an electron 
on the first level of the augmented space in the potential 
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J7i — S. Taking into account the explicit form of C/i ^ S 
with respect to spin- up and spin-down subspaces, 



71 - Ji5_ 
JiS+ 71 - 



where operators 71 and Ji are 

71 = ^i{\s){d\ + \d){s\}, 
Ji = Ji\s){s\, 



one gets 






Gl = 




1 - 


g{ = 




1 - 




-- gIj, 



(B5) 



(B6) 



where 



9l''^ = 



and n+ — S-S+, fi- = S+S-. The physical meaning of 
these formulae is rather transparent. The Green's func- 
tion 5i corresponds to the propagation of the electron 
in a spin conserving part of the potential [/i — E which is 
71 — E^'^), while gY^^ corresponds to scattering on the 

potential 71 — E^'-'--' + Jin±g{^^"' Ji, renormalized with 
respect to the spin-conserving potential due to the inter- 
action with surface magnons on the interface. 

Coming back to the evaluation of scattering matrix 
element toQ, let introduce the "denominator" D, corre- 
sponding to the whole augmented potential [/ — E 



D = 



1-G U- 



(u-t) 



Dqo Dqi 
Dio Dii 



Again, using the technique of the inversion of a matrix 
in the block form and taking into account the elements 
U (|B3|) with respect to the auxiliary space <&, the blocks 
of D can be expressed in terms of the propagator Gi as 
follows 

^00 - [i-g(c/o-E + AGiA)]"\ 

Dio = GiADoo, Doi = Go^Du (B7) 

Dii - fl + GoAGiA) i^oo, 



where we define the propagator Gq, 



Go = 



1-G 



E-^ AGiA 



G = 



Gq Gq 



corresponding to the propagation of the electron in the 
effective potential W —Yj, where 



l^ = [/o + AGiA 




(B8) 



The potential W can be regarded as a renormalization of 
the " virtual" crystal potential C/q of the zero level of the 
augmented space, representing the average of the random 
potential on a site. The renormalization comes from the 
"interaction" A in the auxiliary space with the first level, 
being described by the propagator Gi . Using the explicit 
form of A, 




one can write down the elements and J±: 



w 



+ SGi'^^Uin±gl'^^^ + Sh±G{'"6, (B9) 



U(T)S 



J± 



-/<T(i) ? -i(T)- 
SGi^'^ 



<5Gi^^^Jin±5j^^^(5 



7- 



(BIO) 



One can regard i&Td) as the effective spin-conserving po- 
tential, taking into account the effects of a disorder, and 
J± — as the renormalized electro n-m agnon interaction. 
Finally, the simple algebra using (B3) (B7), yields 



too 




G - E 



Thus, we have obtained the logical result, that the 
t-matrix, averaged over configurations, corresponds to 
scattering on the effective potential W — T,. 

To evaluate the spin-conserving ij)^^^ and spin-flip 
parts of scattering matrix, we introduce the propagators 



m) ^ _ G^ii) (^^ni) _ f;T(i))] GT(i)^ (Bii) 



5o 



corresponding to the potentials w^^ — E^^-^ and define the 
"denominator" 



1-G 



(w -t'j 



Do 



related with the effective potential — E. The elements 
of Do and propagator Go can be expressed via using 
the form of potential W namely 



(j-O 



GT (w^ 



E^ + J+n+g^J_ 
E^ + J^h^g^Jj^ 



G\ 



G^ 



Gq — GlJ+S-gQ, Gq — GgJ-S+gl, 



(B12) 
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and 



1 - [w-^ - + J-h-glj+^ 
gIj+S-4, Do = G^oJ-S+dl, (B13) 
(i) 



where "denominators" dp given by 



"0 — 



Now, using (B12) and (B13) and taking into account the 



obvious relation ioo = {W — S)L'o, one finds 



w' 



According to the ASF, gJ^-*--* and ij*^'' represent the 
configurationaUy averaged quantities, and after the av- 
eraging over magnon degrees of freedom one must get 
{Gl^^^)b = GT(i) and = 0. It means that due to 



the CPA self-consistency conditions the averaged scatter- 
ing t-matrix must vanish and averaged electron propaga- 
tor has to be equal to the effective Green's function ([2l]), 
determined by the self-energy operator E. Carrying out 
the averaging procedure over magnon degrees of freedom 
the similar way as it was proposed earlier in Sec. IIC, it is 
possible to satisfy both of these condition if one assumes 
that 



= LjT(n) + J+(n)n5^(n)J_(n), 
Ti^ = UJ-^{n) + J^{n)ngl{n)J+{n). 



(B15) 



We have pointed out the explicit dependence of matri- 
ces ili^^^'>{n), J±(n) and gl^^\n) on the average magnon 
number n (Eq. p3|) . This dependence is assumed to be 
the same as it comes from the initial definition of opera- 

as functions on n± [Eqs. ([B9|)- 



tors u}'^^^\ J±, and gj^'''^ 



(Bll)], that is the consequence of the adopted approxi- 
mate averaging procedure. The matrices w^*-^-*, J± and 
^q'^' functionally depend on T,^^^\ Due to that, the sys- 
tem (B15) represents the alternative to Eq. (^) form of 
the CPA conditions and it can be simply solved by means 
of successive numerical iterations. 
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